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We show that the dilute Fermi gas quantum critical universality class quantitatively describes the 
Mott /metal crossover of the two-dimensional Hubbard model for temperatures somewhat less than 
(roughly half) the tunneling but much greater than (roughly twice) the superexchange energy. We 
calculate the observables expected to be universal near the transition — density and compressibility 
— with numerically exact determinantal quantum Monte Carlo. We find they are universal functions 
of the chemical potential. Despite arising from the strongly correlated regime of the Hubbard model, 
these functions are given by the weakly interacting, dilute Fermi gas model. These observables and 
their derivatives are the only expected universal static observables of this universality class, which 
we also confirm by verifying there is no scaling collapse of the kinetic energy, fraction of doubly 
occupied sites, and nearest neighbor spin correlations. Our work resolves the universality class of 
the intermediate temperature Mott/metal crossover, which had alternatively been proposed to be 
described by more exotic theories. However, in the presence of a Zeeman magnetic field, we find that 
interplay of spin with itinerant charge can lead to physics beyond the dilute Fermi gas universality 
class. 



Introduction. — The crossover between a Mott insula- 
tor (MI) and a metal in strongly correlated materials is 
an incredibly rich problem pQ. At least two key factors 
play a role: the emergence of a charge gap and the onset 
of magnetic order due to superexchange, and the role of 
each was debated originally by Slater and Mott in the 
1940's and 50's. These intertwined effects — the former 
of which is in general indescribable by a local order pa- 
rameter — give small energy scales at the transition and 
excitations proliferate, and consequently impede a full 
theoretical treatment. A quantum critical perspective 
has frequently been applied to the two-dimensional Hub- 
bard model, the simplest canonical model that includes 
MI and metallic states. At temperatures low compared to 
the superexchange temperature there are numerous pro- 
posed scenarios but few definitive conclusions, and, more 
basically, the phases at zero temperature are unknown. 
Moreover, open questions remain for even the simplest 
quantum critical scenarios (e.g., [2]). 

We address the related question at temperatures above 
the superexchange scale, showing that in an appropri- 
ate regime the "dilute Fermi gas" (DFG) universality 
class quantitatively describes the low energy behavior. 
Although the physics at these temperatures is simpler 
than at low temperatures, it has remained uncertain 
whether and in what regime universal physics manifests, 
and even less clear which universality class describes the 
behavior pQ. Although a common scenario of the doping- 
driven Ml/metal crossover at these temperatures involves 
the DFG theory where the effective number of carriers 
vanishes at the (avoided or fictional) zero-temperature 
transition, such behavior has not been verified with ex- 
act calculations pQ. Other possibilities exist, for exam- 
ple a diverging effective mass at the transition. Another 



possibility comes from zero-temperature determinantal 
quantum Monte Carlo, which indicates a dynamic criti- 
cal exponent z — 4, in contrast to z = 2 predicted by the 
DFG and other low temperature theories This has 
been conjectured to persist up to temperatures above the 
superexchange scale PQ. We show through analysis of 
numerically exact determinantal quantum Monte Carlo 
calculations that the DFG provides a quantitative de- 
scription in an appropriate regime. 

We consider the two-dimensional square lattice Hub- 
bard model defined by the Hamiltonian 

H = E f&i* + E [(£V2K(n- 1) - 1** + hs lz \\) 

with f ia and j ia fermionic annihilation and creation op- 
erators for states at site i and spin state a E {t>l}> 

iMa = ft fi a , n i = "if + n H< s « = n i\ ~ n H> tne sum 
j) indicating a sum over nearest neighbor sites i and 
j, t the tunneling, U the interaction energy, /i the chemi- 
cal potential, and h an external Zeeman field. We defined 
the chemical potential such that the center of the half- 
filling ni = 1 MI at t = occurs at /j, = 0. 

We use determinantal quantum Monte Carlo 
(DQMC) [I] to compute spin densities n a = (n^j), 
compressibility k = dn/dfi, kinetic energy 
K = -tJ2(j},Mhfj<?) ( sum over neighbors j of arbi- 
trary site i), doubly occupied site fraction D = (n^mi), 
and the nearest-neighbor spin correlator X = (si Z Sj Z ) 
with i and j nearest neighbors. We show that for h = 0, 
small enough t/U, and T in the window t 2 /U <T <t 
(we set Ub = 1), the DFG universality class quanti- 
tatively describes the /it-dependence of the expected 
universal observables (n and k) near the Ml/metal 
crossover. We also confirm that, as predicted by the 
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DFG theory, other observables are non-universal. In 
contrast, we find that the Zeeman field dependence of 
the observables is not well described by the DFG or any 
z = 2 universality class, and perhaps is described by 
another, unknown universal theory. 

Many of these observables can be measured in cold 
fermionic atoms trapped in optical lattices [3] , which are 
already in the T ~ t regime. Thus, one can apply a 
similar analysis to that carried out here. In the present 
region, one can use the comparison to validate the faithful 
emulation of the Hubbard Hamiltonian. As these exper- 
iments achieve somewhat colder temperatures, they will 
be able to explore more exotic quantum criticality inac- 
cessible to DQMC [BJ. Even at the present temperatures, 
these may yield insight into quantum critical spectra and 
dynamics that are unobtainable numerically. 

Scaling. — For conventional quantum phase transitions 
described by symmetry breaking of a bosonic order pa- 
rameter, observables O take the form [7J 0(g,T,...) = 
O r (g,T,...) + T 1 + rf / z - 1 /(-o)v I / (( 5 _ 5c )/ T i/(^ 9 )) 

where g is a relevant coupling, the ". . ." of the ar- 
gument indicate irrelevant couplings, vq and v g are, 
possibly non-trivial, scaling dimensions, d is the spatial 
dimension, and z is the dynamic critical exponent. The 
function O r is analytic and is a singular, universal 
contribution. The definition of "universality" is that the 
only dependence of \&o on the irrelevant variables is to 
change the dimensionful constant in that converts 
rpi/iiZg) i n to a dimensionless constant. These constants 
for z = 1 and z — 2 theories have the interpretation of 
an effective sound speed and mass for the excitations, 
respectively. We have assumed g is the only relevant 
coupling. Generally there can be more than one relevant 
variable, but frequently there are only one or two. 

For conserved observables O such as density and mag- 
netization and their dependence on variables coupling to 
these quantities, the scaling structure simplifies to [5] 

Ofji, h, T„ . .) = OM h, T,. . .) + T** fc^, ^)(2) 

where x = 1 + d/z for O — n a and x = d/z for O = k. 
Note that here we include two relevant couplings, h and 
fx, and these are the only two relevant static couplings for 
the DFG universality class we consider below. Close to 
critical point one can Taylor expand the regular pieces 
n a r and k t and keep only the constant terms, which 
we denote and k^°\ For the Ml/metal crossover, 
the constant term for the regular piece for density and 
compressibility are intuitively expected to be n'°) = 1 
and k' ) = 0, and the calculations below confirm this. 

This scaling structure also applies to some more gen- 
eral transitions than symmetry breaking bosonic ones, 
including the DFG transition of interest here [7J. This 



theory is defined by the Hamiltonian 
Hdfg = ~ A** ) n k<T + ^2 h * ( nk t - n H) 

k,CT ^ ' k 

+ f XI c k,T c Li c p-q,i c k+q,T ( 3 ) 

P,k,q 

where m* is the effective mass, g the effective interac- 
tion, n* the effective chemical potential, h* the effective 
magnetic field, c^J are fermionic annihilation (creation) 
operators for momentum k and spin a, and p\^ a = c ] { (;cr c kcr . 
For the present case, it can be shown that /i* = — fi c 
and h* = h — h c . The interaction term is the only rel- 
evant coupling in d < 2 and is marginal in d — 2. We 
stress that this is a (candidate) effective model only near 
the quantum critical crossover, and there is no simple 
relationship between the renormalized parameters g and 
m* and the microscopic parameters U, fj,, h, and t. Our 
DQMC calculations are for the Hubbard model, Eq. 

Despite the fairly simple nature of the Hamiltonian 
in Eq. ([3|, it displays a quantum phase transition with 
much of the phenomenology of general quantum critical 
behavior. The T — phase transition is tuned by [i and 
occurs at /i c = 0. For fi* < 0, the zero temperature 
system is a vacuum of no particles, and at temperatures 
T <C fi* consists of a dilute classical gas. For // > 0, 
the zero-temperature system is a Fermi liquid, and low 
finite temperatures T <C |^*| add dilute quasiparticle ex- 
citations. At finite temperatures satisfying T > 
an intervening "quantum critical" region occurs, and the 
spatial separation between excitations is comparable to 
their thermal de Broglie wavelength. 

If we ignore the marginal interaction term, the DFG 
observables are simply those of a non-interacting gas: 

n^,h,T) = g^log(l + e ^+^) 

m* / I 1 V 

' } ~ 2ntf \eP^~ h ) + 1 + eKv+V + 1 / } 

with f3 = 1/T, and a in the exponential takes values 
{+1, —1} for a = We expect the interactions to 

give logarithmic corrections in d = 2, but we find that 
these are negligible in the regime considered herein. 

Because the density is a derivative of the free energy 
with respect to fj,, & relevant coupling, it will have a sin- 
gular, universal contribution. Indeed, the only universal 
static observables in the DFG theory are the density and 
magnetization (and their derivatives with respect to fi 
and h). Other observables, for example the double occu- 
pancy or nearest neighbor spin correlations, are deriva- 
tives of the free energy with respect to irrelevant cou- 
plings and thus have no singular, universal contribution. 
Consequently, these other observables are non-universal. 

Results — We compute several observables, enumerated 
below, with DQMC for a L x L lattice with L = 10, 
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FIG. 1: (Color online) Universal /^-dependence of the density 
and compressibility for several tunneling rates and interaction 
strengths (t/U's). Data are for T = t/3, except for U/t — 12 
for which T = t/2.5. We observe excellent universal collapse 
for all data at U/t = 6, 8, 12 and also for all but very low hole 
densities at U/t = 4. Top: (n - 1)£/T vs. (jj, - fx c )/T for 
h = 0, where £ = m* /(2%h 2 ). Bottom: k£ vs. (n - 1)£/T. 
We omit £'s in plot labels throughout, as is conventional. We 
fit the non-universal scaling factor £ to obtain collapse with 
the dilute Fermi gas critical theory prediction, Eq. Q. 



discretizing imaginary time evolution into Trotter steps 
of size St — l/(12i), running 1000 Monte Carlo equi- 
libration steps, followed by 100, 000 steps for statistical 
sampling, for each value of n and h. The one exception 
was the t/U = 1/12, T/t = 1/3 calculation, for which 
500, 000 statistical sampling steps were necessary to get 
sufficiently accurate results. Typical statistical error bars 
for the density, not shown, are less than a tenth of a 
percent, smaller than the point size. We have checked 
convergence in each of these parameters. We found that 
the results were unchanged roughly within statistical er- 
ror bars going from L = 8 to L = 16 for the t/U = 1/4 
case where we expect the finite-size effects to be most im- 
portant. We similarly found no change within statistical 
error bars from decreasing the Trotter step to l/(16i). 

To determine /i c , n£\ and k^°\ we plot n/T versus /i 
(for the h = calculations). Eqs. |2]) imply that for 
observables in the scaling region, the resulting curves 
should cross at n = (i c with the crossing value giv- 
ing £K°).From this, we find that = 1, «(°) = 0, 
and fi c /t = {0,-0.83,-1.2,-3} for U/t = {4,6,8,12}. 
We also note that plotting k versus n/T gives a univer- 
sal curve, so one can circumvent determining n c . Note 
that the U/t — 4 value of // c is likely non-zero, but 
is zero within our error bars. The fi c /t used in the 



scaling analysis below are slightly different than these, 
namely n c /t = {—0.3,-0.8,-1.5,-3.2}. These are in 
rough agreement with the crossing points, but shifted by 
5(/j c /t) s» 0.2. These are found by adjusting to get best 
agreement with DFG theory. The difference between the 
two methods is likely due to the fact the crossing is in- 
fluenced by data from non-universal temperatures. It is 
worth comparing these values of fi c /t to expectations in 
the small- and l&ige-t/U limits. Note the trend of de- 
creasing fi c for increasing t: as expected, the MI regime 
shrinks as t increases. For small t/U, we can compare 
this to naive expectations from a simple model of Hub- 
bard bands with width 8t that u c /t = —U/(2t) + 4, 
which gives a c /t = —2 for our smallest t/U = 1/12, 
in rough agreement with our observations. For large 
t/U, the RPA treatment of the antiferromagnet gives 
fj, c « -te~ 27: Vt/u. For our smallest t/U = 1/4, this 
is n c /t = —0.04, also in rough agreement with our ob- 
servations, although this data is far from the large t/U 
limit. 
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FIG. 2: (Color online) Universal pi-dependence of the densi- 
ties and compressibility. Universality breaks down outside of 
t 2 /U < T < t window. Left column: U/t = 6, Right column: 
U/t — 8. Although the collapse for U/t = 6 is fairly good at 
the lowest temperatures, the collapse for U/t = 8 is better for 
a wider range of temperatures, as expected since the universal 
temperature window is larger. 

Figure [T] shows that the density and compressibility 
are universal functions of n for temperatures T = t/3 
(all but t/U = 1/12) or t/2.5 (t/U = 1/12, because 
it is difficult to get accurate compressibilities at lower 
temperatures). Each panel shows shows the appropri- 
ately rescaled observables for t/U = 1/4,1/6,1/8, and 
1/12. The collapse indicates z = 2 scaling behavior for 
several t/U's at t/T = 1/3. Furthermore, Fig. [I] shows 
quite remarkably that the non-interacting DFG expres- 
sions, Eq. Q, quantitatively describe the scaling func- 
tions of the Hubbard model in the strongly correlated 
regime near the Ml/metal crossover. Although in d = 2 
we expect logarithmic corrections to the non-interacting 
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scaling functions, these are unobservably small. 

Fig. [2] gives further evidence for this by showing col- 
lapse of data at multiple temperatures for t/U = 1/8 
and t/U = 1/6, and we find that the DFG quantitatively 
describes the //-dependence for temperatures satisfying 
t 2 /U <T<t,a fairly natural result. The T = t/3 
results in Fig. [T] are also in this temperature window. 
We also have checked that, in contrast to the plots of 
rescaled variables, e.g. (n — 1)/T vs (fx — Hc)/T, plotting 
variables without rescaling or plotting rescaled variables 
corresponding to other z's, e.g. z — 1 or z — 4, does not 
give good collapse (see Supplementary Information). 

The data collapse observed for t 2 /U < T < t breaks 
down for both lower and higher temperatures. At higher 
temperatures this is natural because the excitations be- 
gin to probe energy scales beyond the bottom of the ex- 
citation band, where the dispersion is well described by 
6fc = h 2 k 2 /(2m*), and begin to see the microscopic band 
structure. At lower temperatures, the spins are no longer 
independent, and develop magnetic correlations. These 
correlations at high temperatures depend on microscopic 
details governing superexchange physics. Only at much 
lower temperatures, T -C t 2 /U will universal behavior re- 
cover, and it will likely be more exotic, for example one 
of the scenarios referenced in the introduction. 

To compare scaling functions obtained with DQMC 
to the DFG theory's, one requires a single non-universal 
scaling factor, the DFG effective mass in Eq. Q. 
From this, we find the ratios of m*'s for t/U = 
1/4,1/6,1/8,1/12 are 1 : 1.7 : 2.4 : 4.1. The m* in- 
creases as t/U decreases, roughly as 1/t for small t. 

Figures [l] and [2] illustrate that universal scaling is best 
for small values of t/U, becoming quantitative over the 
full visible range only for t/U < 1/8. This is natural, 
since only for small t/U is there a large separation be- 
tween the t 2 /U and t temperature scales that must hold 
in order for the temperature to satisfy t 2 /U « T « t. 

In the DFG universality class, n-j-, n\., and k and 
their derivatives are the only universal static observ- 
ables. Figs. [T] and [2] confirmed their universality. We also 
confirmed that other observables are non-universal: the 
nearest neighbor spin correlations, doubly occupied sites, 
and kinetic energy (see Supplementary Information). 

In contrast to naive expectations, we find that the 
DFG fails to capture the /i-dependence of the spin den- 
sities. The reason is that in addition to the itiner- 
ant charge carriers governed by the DFG, there is a 
Mott background with spin degrees of freedom. For 
t/U -C 1, this background is described by a high tem- 
perature T 3> t 2 /U Heisenberg model, which reduces to 
a single site problem. The spin densities from this back- 
ground are rvf = 1/(1 + e 2 ^ h ) and — 1/(1 + e~ 2 P h ). 
At T = 0, these are singular (step functions) at h = 0, 
indicating a non-analytic contribution to the spin densi- 
ties that is not captured by the dilute Fermi gas theory. 
Thus at fx = jj, c for h = 0, there are coinciding singu- 



larities from this "background" and DFG contributions. 
Rather than simply summing, it appears they may com- 
bine to give new universal behavior, with z w 4. This is 
elaborated in the Supplementary Information. 

Summary. — We have computed the density, compress- 
ibility, fraction of doubly occupied sites, kinetic en- 
ergy, and nearest neighbor spin correlations of the two- 
dimensional square lattice Hubbard model with deter- 
minantal quantum Monte Carlo near the Mott/metal 
crossover. Resolving a long discussion in the literature, 
we find that the behavior is quantitatively described by 
the dilute Fermi gas universality class. For the observ- 
ables that this class predicts to be universal — spin densi- 
ties and compressibilities — we find the dilute Fermi gas 
theory quantitatively gives the scaling functions at h = 
for a range of t/U < 1/6 and temperatures t 2 /U <T <t. 
We also confirm that there is no universality for observ- 
ables the dilute Fermi gas theory predicts to be non- 
universal. For h ^ 0, the singular contribution from the 
non-itinerant spin degrees of freedom invalidates the di- 
lute Fermi gas theory, and interaction of localized spins 
with itinerant carriers may give new universal physics. 
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COMPARISON OF RESCALING FOR z = 1, z = 2, 

AND z = 4 
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FIG. 3: Rescaled compressibility (n/T d/z - 1 ) versus rescaled 
chemical potential (/i — /i c )/T for U/t = 8 and dynamic expo- 
nent z — 1,2,4, top to bottom. The z — 2 rescaling obtains 
good collapse while the 2 = 1 and z — 4 both fail to give any 
significant collapse. The other values of U/t we considered 
show similar results. 



were no good crossings (fj, c = and /i c = —1.5, respec- 
tively, were the best or at least as good as any others). 
Note that in the main text, the presented data used fj, c 
determined to give good collapse of the low-temperature 
data to the dilute Fermi gas theory, and the resulting 
value of fj, c was slightly different than with the crossing 
method (e.g., /i c = —1.5 instead of fi c = —1.2 for the 
z = 2 rescaling of the U/t = 8 data) . Since we have no 
natural theory to use for comparison in the z = 1 and 
z — 4 cases, determining fi c by the crossing points was 
more fair for the comparison of Fig. [3j 



FAILURE OF NON-UNIVERSAL OBSERVABLES 
TO RESCALE 



As discussed in the main text, besides the /x- and h- 
derivatives of the free energy, for example the densities 
and compressibilities, we expect no other static universal 
observables. We have confirmed that other observables 
indeed do not collapse: the kinetic energy K, nearest 
neighbor spin correlator X, and double occupation prob- 
ability D. This remains the case, even if we allow the 
non-universal constant, D(T — 0, fi — fi c ), and the scal- 
ing exponent to take arbitrary values. Fig. [4] illustrates 
the best collapse we could obtain for the kinetic energy 
at t/U = 1/8, found by matching the small hole den- 
sity behavior of the two lowest temperature curves. As 
expected, the collapse is very poor. These results are 
typical. 
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Figure [3] illustrates that the z — 2 rescalings offer a 
good universal collapse of the data, while other candi- 
date z's fail. We show rescalings associated with z = 1,2, 
and 4, and only the z = 2 case provides an adequate de- 
scription of the data. We determined \i c from the best 
crossing point of n/T d l z , but for z = 1 and 2 = 4, there 



FIG. 4: Rescaled kinetic energy (K + 0.87)/r l ei versus 
rescaled density (n — 1)/T, for t/U = 1/8. We have adjusted 
the offset (0.87) and the exponent (1.6) to find the best col- 
lapse for the two lowest temperature curves, but collapse is 
still poor. 
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^-DEPENDENCE OF OBSERVABLES 
DISAGREES WITH DFG THEORY 

We have noted that \i- and fr-are (the only) relevant 
tuning parameters within the dilute Fermi gas (DFG) 
theory. Thus we might expect the DFG to govern the 
universal /i-dependence of observables. However, we 
have computed the spin-up and spin-down densities of 
the Hubbard model as a function of h for // = /i c and 
t/U = 1/8 and find that the DFG theory does not ad- 
equately describe these observables. More strongly, no 
scaling collapse is obtained for z — 2, ruling out not only 
the DFG universality class, but any z = 2 theory. Fig. [5] 
(top panel) demonstrates the lack of z = 2 scaling, and 
the middle panel shows the DFG theory curves for com- 
parison. 

To understand this lack of collapse, we reconsider the 
physical picture leading to the application of the dilute 
Fermi gas theory. We also note that no formal justifi- 
cation has to this point been available, so that we must 
rely on our physical intuition. The basic picture was 
of itinerant dopants governed by the DFG theory cou- 
pled to an essentially inert Mott background. The Mott 
background's only effect was to renormalize the DFG pa- 
rameters and provide a non-universal, non-singular con- 
tribution to the observables, which was irrelevant near 
the quantum critical point. 

However, when a magnetic field is applied, treating 
the Mott background as inert is dangerous. To see this, 
consider the case where there are no itinerant charge car- 
riers and thus the DFG plays no role. The Mott back- 
ground may be described by a Heisenberg model, and 
since we are in the high temperature limit T ^> t 2 /U, 
this reduces to a single site problem. In this case, the 
spin densities are = er^ h / {er^ h + e' 3 ' 1 ) and nj, = 
e@ h I {e~P h + e^ h ). We see these are singular (step func- 
tions) at h = and T = 0, so that the "Mott back- 
ground" gives a singular contribution here, rather than a 
regular contribution as in the h = case, tuning ji. This 
indicates that the singular contribution will include both 
a universal contribution from the DFG theory and from 
the "background." The extent to and way in which these 
two pieces of physics couple at h ^ is a complicated 
issue that we leave to future work. 

Interestingly, we do observe rather good collapse for 
z = 4, reminiscent of the z = 4 exponent found at T = 
where magnetic physics, including superexchange, is cer- 
tainly relevant 1 . However, we emphasize we are well 
above the superexchange temperature. It would be in- 
triguing to understand this result better, and in partic- 
ular if it indicates a new universality class governing the 



t 2 jU <C T -C t Hubbard model. Presumably such a uni- 
versality class would consist of free spins coupled to a 
DFG in some way, although to be consistent with our 
results in the main text it must reduce to the DFG uni- 
versality class when h = 0. 
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FIG. 5: Top: determinantal quantum Monte Carlo data for 
spin-down density (top curves) and spin-up density (bottom 
curves) rescaled at U/t = 8, assuming a dynamical critical 
exponent z = 2. Going from "inner curves" to "outer curves," 
the temperatures are temperatures T/t — 2,1,1/2,1/3. No 
scaling collapse is observed. Middle: dilute Fermi gas theory 
predictions for comparison. Bottom: same as top, but results 
rescaled according to z = 4. 
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